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1 Introduction 

Deep imaging in the rest-frame UV has constrained both the evolution of the 
cosmic star formation rate density [1] and its time integral, the growth of the cos- 
mic stellar mass density [2] . Short-wavelength studies give an incomplete picture, 
however, since an important population of high-redshift galaxies is heavily dust- 
obscured. The strength of the extragalactic mid- and far-IR/submillimeter back- 
ground indicates that about half of the cosmic energy density comes from dusty 
luminous and ultra-luminous infrared galaxies (LIRGs/ULIRGs: Lir, ~ 10 115 to 
10 13 ' 5 Lq) at z > 1 [3,4]. Because the brightest of these submillimeter galaxies 
(SMGs; see [5] and references therein) tend to lack strong X-ray emission [6], 
their large IR luminosities probably correspond to high star formation rates [7] . 
As the strikingly different appearances of the Hubble Deep Field at 0.83 /im [8] 
and 850 /im [9] exemplify, SMGs are rarer and forming stars much more intensely 
than typical optically selected systems [8,9]. Here we discuss new observations 
that shed light on the importance of SMGs in the history of galaxy mass assembly 
(all numbers assuming a flat fl^ = 0.7 cosmology with H = 70kms~ 1 Mpc -1 ). 



2 Observations: PdBI millimeter interferometry 



Since SMGs in most cases have relatively poorly known positions and frequently 
have only weak counterparts in the rest-frame UV and optical, redshifts and 
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PdBI CO (3-2/43) Survey 
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Fig. 1. Integrated CO spectra of the seven detected PdBI survey sources [11,12] plus 
SMMJ 02399-0136 [13]. Source names and redshifts are listed in the insets; flux density 
and velocity scales are the same for all sources. 



spectroscopic parameters were initially confirmed with CO interferometry for 
only two of the few hundred detected sources [5]. Recently a subset of the 
authors have initiated a new program to obtain spectroscopic redshifts with 
the Keck telescope for a number of sources detected with SCUBA at 850 /im, 
in most cases aided by more precise positions derived from deep 1.4 GHz VLA 
maps [10]. Here we report the first results of the millimeter follow-up of these 
SMGs, undertaken with the IRAM Plateau de Bure Interferometer (PdBI). In 
the B, C, and D configurations of the PdBI, we searched for the CO(3-2) and 
CO (4-3) rotational transitions (redshifted to the 3 mm atmospheric window) in 
a sample of so far twelve SMGs with 5*850 > 5mJy and reliable UV/optical 
spectroscopic redshifts. Integration times per source ranged from 13 to 30 hours, 
for a total integration time of about lMs in the 2002/2003 observing season. We 
have detected significant CO emission in 7 out of the 12 sources, thus roughly 
quadrupling the number of CO-confirmed SMG redshifts [11,12]. Figure 1 shows 
the integrated CO(3-2/4-3) spectra for the seven new sources and for the z — 
2.81 system SMM J02399-0136 [13]. Figure 2 shows the integrated CO maps 
superposed on rest-frame UV (/-band) and optical (if-band) images of three of 
the sources [11]. 

Our PdBI observations confirm the median z ~ 2.4 determined from rest- 
frame UV/optical spectroscopy [10], although in several cases, the UV red- 
shift is inferred from a source that is physically distinct from the SMG (e.g., 
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Fig. 2. Integrated CO maps of three SMGs studied with the PdBI [11], superposed on 
greyscale images of rest-frame UV (center panel) and optical (left and right panels) 
emission. FWHM synthesized beams are shown as hatched ellipses. In the left panel, 
the asterisk marks the position of the ERO N4 (uncertainty ±0.5" - see text), while 
the edge-on spiral galaxy Nl (2" SE of the CO source) lies in the foreground cluster 
at z = 0.18. In the center and right panels, filled squares mark millimeter continuum 
positions. Arrows show the positions of radio sources H6 and H7; the stronger of these 
(H6) is a narrow-line Seyfert galaxy with Lya indicating z = 3.349. CO(4-3) emission 
from N2 850.4 remains largely unresolved. 



SMM J09431+4700 in Fig. 2). SMGs likely sample very dense environments that 
may harbor additional, UV-bright star-forming galaxies. 

SMGs are rich in molecular gas. Correcting for foreground lensing wherever 
appropriate (and known), and using the M gas / Lqo conversion factor appropriate 
for z ~ 0.1 ULIRGs, the median gas mass of the eight SMGs in Fig. 1 (plus the 
z = 2.56 galaxy SMM J14011+0252 [14,15]) is 2.2 x 10 10 M Q . This value is 
similar to the median molecular gas mass found in high-z QSOs and about three 
times greater than in local ULIRGs. The inferred median far-IR luminosity of 
the sample (again corrected for lensing) is about is 1.3 x 10 13 L©, corresponding 
to a star formation rate of about 1100 M©yr -1 . 

Our most important finding is that the brightest SMGs (Ssso > 5mJy), 
contributing about 25% of the submillimeter background, have large dynami- 
cal masses. The median velocity dispersion of our sample is 255(±38) kms -1 . 
The ratio of the squares of the median velocity dispersions then implies that 
our SMGs have dynamical masses at least 13(±4) times greater than the z ~ 3 
Lyman break galaxy (LBG) population [16], and at least 5(±1.5) times greater 
than optically selected populations at z ~ 2 ([17] - see also this volume). There 
are indications that SMGs are larger than LBGs in which case these estimates 
are lower limits [13,11,18]. Assuming a rotating disk geometry at an inclina- 
tion of 40°, the median velocity dispersion corresponds to a dynamical mass of 
1.1 x 10 11 Mq within R — 4kpc. Based on its spatially resolved rotation curve, 
SMM J02399-0136 has a total dynamical mass of 3 x 10 11 M Q within R = 8kpc, 
half of which is probably stellar [13]. Near-IR photometry indicates that the me- 
dian stellar mass for the z ~ 3 LBG population is 1.1 x 10 10 .Mq[19], consistent 
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with the dynamical mass ratio estimated above. Dynamical and stellar masses 
of the SMGs in our sample are thus comparable to those of local m* galaxies. 

3 Observations: SPIFFI NIR integral field spectroscopy 
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Fig. 3. J, H, and K spectra of J14011 extracted from a 2" x f.5" aperture [22]. 
Superposed is the model STARS spectrum for a 200 Myr old continuous star formation 
episode with solar metallicity, a f — 100 M Q Salpeter IMF, and extinction Av = 0.7. 

We have also conducted a detailed rest-frame optical case study of the bright 
source SMM J1401 1+0252 (hereafter J14011), a z = 2.565 SMG lying behind the 
z = 0.25 cluster A1835 [20]. For this study we used SPIFFI, the SPectrometer 
for Infrared Faint Field Imaging [21], at the ESO VLT. In three runs between 
February and April 2003, we observed J14011 in the J, H, and K bands, with 
a FWHM spatial resolution of about 0.6" and spectral resolutions ranging from 
125 to 200 km s -1 . SPIFFI obtains simultaneous spectra for each of 32 x 32 
contiguous 0.25" pixels. The total on-source integration times were 60 minutes 
in J, 95 minutes in H, and 340 minutes in K [22]. 

Within the bright eastern Jl component of J14011, the SPIFFI data reveal 
substantial spatial variation in Ha equivalent width: nebular line emission is 
extended over about 1" (lensing corrected) and is strongest away from the central 
continuum peak. The integrated spectrum of Jl (Fig. 3) appears to exhibit a 
continuum break between the observed J and H bands. The feature's most likely 
origin is the Balmer break of a > 100 Myr old stellar population at z ~ 2.565. The 
Jl line fluxes allow us to measure the R23 estimator for its oxygen abundance. 
After correcting for the E(B — V) = 0.18 implied by the Balmer decrement, we 
derive an oxygen abundance of 8.96 ±0.10, or 1.9 times solar. The observed high 
[Nil] /Ha ratio implies that J14011 cannot fall on the lower-metallicity branch of 
the [O/H] vs. i?23 relation. In contrast to the high abundances claimed for some 
high-z QSOs, this metallicity refers to the entire galaxy on ~ lOkpc scales. 

For a closed-box model with solar yield, this metallicity corresponds to a 
gas/baryonic mass fraction of cxp (— Z/Zq) ~ 0.16. From its intrinsic (corrected 
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for a lensing magnification of 5) gas mass of ~ 1.3 x 1O 1O M0, we therefore 
infer a total baryonic mass ~ 8 x 10 10 M© and stellar mass ~ 7 x 1O 1O M0. 
This last value is empirically supported by the location of J 140 11 in the local 
mass-metallicity relation. Adopting a recent version of that relation [23], we 
would expect a stellar mass of ~ 6 x 10 10 M©. From the closed-box model with 
continuous star formation at the current rate (380 M© yr _1 ) we infer the overall 
duration of the star forming activity in J14011 to be tsF — 220 Myr. This age 
is entirely consistent with the strength of the Balmer break seen in Fig. 3. As 
a further consistency check, we can estimate a dynamical mass for J14011 from 
the ~ 2.2" and ~ 180 km s -1 separations between Jl and [second component] 
J2, assuming these orbit around a common center. For a predominantly north- 
south lensing shear, the result is ~ 3.4 x 10 10 sin~ 2 i M©, consistent with the 
moderately low inclination also implied by the Ha morphology. 

4 SMGs and mass assembly 
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Fig. 4. Comoving number densities of galaxies with baryonic masses > 10 11 Mq as a 
function of redshift [22]. Triangle and open squares show densities of massive stellar 
systems at z = [25] and z ~ 1 [26]; circle shows density for massive SMGs at z ~ 2.7, 
with a factor of 7 correction for burst lifetime. Solid curves show the predictions of semi- 
analytic modelling by the "Munich" and "Durham" groups [27,28]; dashed curves show 
the corresponding number densities of halos with available baryonic masses > 10 11 Mq. 
The two models use the same halo simulations but assume different J?b. 

Our IRAM and VLT results on mctallicity and gas/dynamical mass build 
a compelling case that the most luminous SMGs (Sg5o ~ 5 to lOmJy) have 
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dynamical and - within the factor 2 uncertainties of the present estimates - also 
stellar masses of ~ 10 11 M Q . This stellar component has plausibly formed in 
the current starburst of duration ~ 100 Myr. How does the inferred comoving 
volume density of SMGs compare to theoretical models of galaxy formation? 

From the observed surface density of SMGs with S$5o > 5 mJy [24] , we 
infer that radio-detected SMGs with baryonic masses > 10 11 Mq in the range 
1.8 < z < 3.5 have a comoving number density of 1.6tl'% x lCT 5 Mpc~ 3 [22]. 
This density needs to be corrected upward to account for sources that are in 
the same mass bin but no longer luminous enough to have been detected at 
850 (im. Relative to the 1.8 Gyr elapsed over the range 1.8 < z < 3.5, the ~ 
220 Myr age and ~ 30 Myr gas exhaustion timescale for the starburst in J14011 
imply a lifetime correction factor of about 7. Figure 4 plots this prediction, 
together with the comoving number densities of galaxies with stellar masses 
above the 10 11 M© threshold at z ~ [25] and z ~ 1 [26]. Also plotted are 
the predictions of two semi-analytic models [27,28]. Since all theoretical and 
observational values assume the same cosmology and a consistent (either Miller- 
Scalo or 1 — 100 M© Salpeter) IMF, we are secure in concluding that the models 
underpredict the number densities of massive galaxies at z ~ 2.5 [22]. Couching 
Fig. 4 in terms of mass emphasizes that the observed surface densities of SMGs 
cannot be explained using a very flat IMF alone. We also note that the models 
contain enough dark halos of total mass > 10 11 (Qm/ M Q to account for 
the observed comoving number densities of massive SMGs, provided all of their 
baryons can be rapidly and efficiently assembled into galaxies. SMGs thus impose 
a powerful "baryonic mass assembly" test at the upper end of the galaxy mass 
function [13]. 

Our correction of the observed number density of SMGs to account for their 
less dust-luminous descendants raises the question of what types of galaxies 
SMGs can plausibly evolve into. Although LBGs are clearly excluded due to 
their smaller masses, an intriguing alternative is the newly identified population 
whose red J — K colors can stem from a strong Balmer break at z ~ 2.5 [29]. 
These galaxies appear to be strongly clustered [30] and - assuming a uniform 
distribution in the redshift range 2 < z < 3.5 [31] - have a comoving number 
density ~ 1.8 x 10~ 4 Mpc~ 3 . Relatively few sources with red J — K colors are 
also luminous (sub)millimctcr emitters, indicating that the two populations have 
little direct overlap. With three examples in the HDF-S having (for a Miller-Scalo 
IMF) stellar masses (0.6 - 1.4) x 10 11 M© [32], it would seem plausible that the 
more luminous objects with red J — K colors at this epoch could have passed 
through a SMG phase at higher redshift. 
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